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\Q , Abstract 

We present detailed time-averaged X-ray spectroscopy in the 0.5-10 keV band of the Seyfert 1.9 galaxy 
■ NGC 2992 with the Suzaku X-ray Imaging Spectrometers (XIS). The source had a factor ~ 3 higher 2- 

Ph| 10 keV flux (~ 1.2 x 10 erg cm s ) than the historical minimum and a factor ~ 7 less than the 

q ■ historical maximum. The XIS spectrum of NGC 2992 can be described by several components. There is 

S-^ | a primary continuum, modeled as a power-law with a photon index of T = 1.57i "o3 that is obscured by a 

Compton-thin absorber with a column density of 8.0^Q5 x 10 21 cm" 2 . There is another, weaker, unabsorbed 
power-law component (modeled with the same slope as the primary) , that is likely to be due to the primary 
continuum being electron-scattered into our line-of-sight by a region extended on a scale of hundreds of 
parsecs. We measure the Thomson depth of the scattering zone to be r es = (0.073 ± 0.021)/[fi/47r], where 
fl/An is the fraction of the sky covered by the zone (as seen from the X-ray source) that is visible to the 
observer. An optically-thin thermal emission component, which probably originates in the same extended 
region, is included in the model and yields a temperature and luminosity of kT = 0.656^0 06i keV and 
~ 1.2 ± 0.4 x 10 40 erg s _1 respectively. We detect an Fe K emission complex which we model with broad 
and narrow lines and we show that the intensities of the two components are decoupled at a confidence 
level > 3<t. The broad Fe Ka line has an equivalent width of 118j^g 2 eV and could originate in an accretion 
disk (with inclination angle greater than ~ 30°) around the putative central black hole. The narrow 
Fe Ka line has an equivalent width of 163^26 °V an d is unresolved (FWHM < 4090 km s" 1 ) and likely 
originates in distant matter. The absolute flux in the narrow line implies that the column density out of 
the line-of-sight could be much higher than measured in the line-of-sight, and that the mean (historically- 
averaged) continuum luminosity responsible for forming the line could be a factor of several higher than 
that measured from the data. We also detect the Fe K/3 line (corresponding to the narrow Fe Ka line) 
with a high signal-to-noise ratio and describe a new robust method to constrain the ionization state of Fe 
responsible for the Fe Ka and Fe K/3 lines that does not require any knowledge of possible gravitational 
and Doppler energy shifts affecting the line energies. For the distant line-emitting matter (e.g. the 
putative obscuring torus) we deduce that the predominant ionization state is lower than Fe vin (at 99% 
confidence), conservatively taking into account residual calibration uncertainties in the XIS energy scale 
and theoretical and experimental uncertainties in the Fe K fluorescent line energies. From the limits on a 
possible Compton-reflection continuum it is likely that the narrow Fe Ka and Fe K/3 lines originate in a 
Compton-thin structure. 
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1. Introduction 

The Fe K emission line profiles in active galactic nu- 
clei (AGN) generally consist of a prominent core centered 
at - 6.4 keV (typically with FWHM < 10,000 km s" 1 or 
so) and in some sources an additional underlying broader 
component. The relatively narrow core at ~ 6.4 keV 
is found to be ubiquitous in nearby, low to moderate 
luminosity AGN (e.g. Nandra et al. 1997a, b; Sulentic 
et al. 1998; Weaver, Gclbord, & Yaqoob 2001; Reeves 
2003; Page et al. 2004; Yaqoob & Padmanabhan 2004; 
Jimcncz-Bailon et al. 2005; Zhou & Wang 2005; Jiang, 
Wang, & Wang 2006; Lcvenson et al. 2006). The broad 
Fe K line component potentially carries signatures of the 
inner accretion disk in the vicinity of the putative black 
hole by virtue of gravitational and Dopplcr energy shifts 
(e.g. Tanaka et al. 1995; Nandra et al. 1997a, b; Reynolds 
1997; Turner et al. 1998; Iwasawa et al. 1999; Yaqoob 
et al. 2002; Reynolds & Nowak 2003; Fabian & Miniutti 
2005). The bulk of the Fe K line core is likely to originate 
in matter located much further from the central black hole 
than the outer disk (e.g. the outer parts of the optical 
broad line region, the putative parsec-scale torus struc- 
ture, or the optical narrow line region). In addition to the 
two principal components, there are an increasing number 
of reports of cases of narrow lines at energies higher than 
6.4 keV (from ionized and/or outflowing Fe - see Bianchi 
et al. 2005, and references therein) and lower than 6.4 keV 
(from apparently rcdshiftcd localized emission - Turner 
et al. 2002; Iwasawa, Miniutti, & Fabian 2004; Turner 
et al. 2006, and references therein). 

Deconvolving the disk and distant-matter components 
of the Fe K line is extremely challenging because the disk 
can also contribute to the line core emission (from radii 
greater than ~ 100r g = 100GM/c 2 or so). In other words, 
line emission from disks observed at low inclination angles 
(relative to the disk normal) is partially degenerate with 
line emission from distant matter. Line emission in the 
wings of the profile from regions close to the black hole is 
smeared out over energy and therefore difficult to detect 
against the continuum. In order to derive the most reli- 
able parameters and constraints for the disk and distant- 
matter Fe K line components, one must simultaneously 
model both, otherwise the model parameters will be biased 
and misleading (e.g. see discussions in Weaver & Reynolds 
1998, Yaqoob & Padmanabhan 2004). So far only some of 
the brightest AGN have sufficient signal-to-noise ratio to 
produce useful constraints from dual-line models (e.g. see 
Yaqoob et al. 2002, Yaqoob & Padmanabhan 2004, and 
references therein) . For the brightest sources the Fe K line 
core is usually unresolved even with the Chandra gratings 
because high spectral resolution and high throughput are 
required. Therefore, the location of the "distant" Fe K 
line emitter remains unknown. Nevertheless, decoupling 
the intensities of the broad and narrow lines (in the sense 



of showing that the data require both to be non-zero at 
some confidence level) would be an important step in con- 
straining models of the central engine in AGN. Such de- 
coupling has only been possible for a few of the brightest 
AGN (e.g. see Fabian & Miniutti 2005). We also note 
that detection of the broad Fe K line in type 1.5-2 AGN 
is even more difficult than in type 1 objects because the 
X-ray spectra of the former class are more complex and 
are usually affected by significant intrinsic X-ray absorp- 
tion (e.g. see Turner et al. 1998; Weaver & Reynolds 1998; 
Matt et al. 2000; Guainazzi et al. 2001; Risaliti 2002). 

Almost all of our knowledge of the Fe K line emission 
in AGN comes from the Fe Ka line. Although Fe K/3 
line emission has been reported in a few type 1.5-2 AGN 
(Sambruna et al. 2001; Molendi, Bianchi, & Matt 2003; 
Matt et al. 2004), constraints on the line energy and equiv- 
alent width (EW) have been loose because the branching 
ratio, Fe K/3/Fe Ka, is less than 0.145 (Palmeri et al. 2003, 
and references therein). Yet the Fe K/3 line carries addi- 
tional information on the ionization state of line-emitting 
Fc that would supplement the information carried by the 
Fe Ka line. 

In this paper we present the results of detailed X-ray 
CCD spectroscopy of the nearby (z = 0.00771, Keel 1996) 
Seyfert 1.9 galaxy NGC 2992 using data from Suzaku ob- 
servations. We show that the intensities of the broad and 
narrow Fe Ka lines are decoupled for the first time in this 
source. We also detected the Fe K/3 line corresponding 
to the Fe Ka line core and present measurements of its 
ccntroid energy, intensity, and EW. We describe a new 
method that utilizes the Fc K/3 line to deduce tight limits 
on the ionization state of Fe in the matter responsible for 
emitting the line core. For the Fe Ka and Fe K/3 line anal- 
ysis we modeled the 0.5-10 keV spectrum self-consistently 
so we also report the broadband results. 

NGC 2992 has been observed by every X-ray astronomy 
mission since the time it was discovered by HEAO-1 to be 
one of the brightest hard X-ray AGN in the sky, with a 2- 
lOkeVfluxof ~ 7.2- 8.6 x lO" 11 erg cm" 2 s" 1 (Piccinotti 
et al. 1982). In more than a quarter of a century of X- 
ray observations, the hard X-ray flux has varied by over a 
factor of 20, corresponding to a range in the intrinsic 2- 
10 keV luminosity of ~ 0.55—11.8 x 10 42 erg s _1 (assuming 
i7 = 70kms^ 1 Mpc" 1 , A=0.73, Q= 1). During an ASCA 
observation in 1994, NGC 2992 was in its lowest contin- 
uum flux state thus far observed, with a 2-10 keV flux of 
~4x 10~ 12 erg cm~ 2 s _1 (Weaver et al. 1996). During 
that observation the Fe K line equivalent width (EW) was 
very high (~ 500 — 700 eV), indicating that the line inten- 
sity had not responded to the declining continuum. In 
1997 a BeppoSAX observation found NGC 2992 to be still 
in a low continuum state, and the Fe K line EW was even 
higher, but then a second BeppoSAX observation in 1998 
revealed that the source continuum had fully recovered to 
its bright HEAO-1 state (see Gilli et al. 2000). In this 
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state, the EW of the Fe K line was very low (< 100 eV) 
since the line intensity again did not respond to the change 
in continuum level. These measurements were mostly sen- 
sitive to the Fe K line core, which from the lack of vari- 
ability, likely arises in a distant, parsec-scale reprocessor, 
as suggested by Weaver et al. (1996). 

The results from an XMM-Newton observation of 
NGC 2992 in May 2003 are unpublished. We examined 
those data and found that they suffered from heavy pile- 
up. A detailed analysis of the XMM-Newton data is be- 
yond the scope of this work and will be reported elsewhere. 

NGC 2992 was the subject of an RXTE monitoring 
campaign consisting of 24 observations in the period 2005 
March to 2006 January. The Suzaku observations that 
are the subject of the present paper were made during the 
RXTE campaign and were quasi-simultaneous with some 
of the RXTE observations. The scope of the present paper 
is restricted to the Suzaku data only, but a detailed analy- 
sis of the data from the RXTE campaign will be presented 
elsewhere. 

The present paper is organized as follows. In section 2 
we describe the observations and details of the data reduc- 
tion. In section 3 we describe the spectral fitting analysis 
and in section 4 we discuss the detailed results of that 
analysis, including brief comparisons with historical data 
and results. In section 5 we describe a new method to 
tightly constrain the ionization state of Fe using the nar- 
row Fe Ka and Fe K/3 emissions lines. In section 6 we 
summarize our results and give our conclusions. 

2. Observations and Data Reduction 

The joint Japan/US X-ray astronomy satellite, Suzaku 
(Mitsuda et al. 2006), was launched on 10 July, 2005. 
NGC 2992 was observed on three occasions in 2005, 
November and December (the observation log is given in 
table 1). Suzaku carries four X-ray Imaging Spectrometers 
(XIS - Koyama et al. 2006) and a collimated Hard X-ray 
Detector (HXD - Takahashi et al. 2006). Each XIS con- 
sists of four CCD detectors at the focal plane of its own 
thin- foil X-ray telescope (XRT - Serlcmitsos et al. 2006), 
and has a ficld-of-vicw (FOV) of 17.8' x 17.8'. One of 
the XIS detectors (XIS1) is back-side illuminated (BI) 
and the other three (XIS0, XIS2, and XIS3) arc front- 
side illuminated (FI). The bandpass of the FI detectors 
is ~ 0.4 - 12 keV and ~ 0.2 - 12 keV for the BI detector. 
The useful bandpass depends on the signal-to-noise ratio 
of the source since the effective area is significantly dimin- 
ished at the extreme ends of the operational bandpasses. 
Although the BI CCD has higher effective area at low en- 
ergies, the background level across the entire bandpass is 
higher compared to the FI CCDs. The spectral resolu- 
tion of the XIS detectors has been continuously degrading 
with time but the current instrument response matrices 
employ the resolution before the onset of degradation, as 
described in Koyama et al. (2006). Although we used 
these standard response matrices for modeling the XIS 
data, we measured the widths of the Mn Ka lines from 
the on-board 55 Fe calibration sources (two per XIS) using 



the actual NGC 2992 observations in order to correctly 
interpret spectral fitting results. Details are given in sec- 
tion 2.3. 

2.1. HXD Data 

The HXD consists of two non-imaging instruments (the 
PIN and GSO - see Takahashi et al. 2006) with a com- 
bined bandpass of ~ 10 - 600 keV. Both of the HXD 
instruments arc background-limited. We found that for 
the NGC 2992 observations the average on-source count 
rate in the full HXD/PIN band (~ 10- 120 keV) was 
~ 0.75 ct/s, whilst the corresponding estimated back- 
ground count rate was 0.65 ct/s so that NGC 2992 con- 
stitutes only ~ 13% of the total signal. The situation 
is even worse for the GSO which has a smaller effective 
area than the PIN. In order to obtain reliable background- 
subtracted spectra, the background spectrum must be 
modeled as a function of energy and time. The back- 
ground model is still under development and cannot yet 
reproduce the background spectral and temporal behav- 
ior with sufficient accuracy to give reliable HXD spectra 
for NGC 2992. Therefore, in this paper we will restrict 
our analysis to the XIS only and defer analysis of the 
simultaneous HXD data to future work when the back- 
ground model systematics are lower. The observations of 
NGC 2992 were in fact optimized for the HXD in terms of 
positioning the source at the aimpoint for the HXD (the 
so-called "HXD-nominal pointing") which gives a some- 
what lower count-rate in the XIS than the "XlS-nominal" 
pointing but gives ~ 10% higher HXD effective area. 

2.2. XIS Data Filtering, Cleaning, and Selection 

The XIS observations of NGC 2992 were made in so- 
called 3x3 and 5x5 edit modes which correspond to 
9 and 25 pixel pulse-height arrays respectively that are 
telemetered to the ground (see Koyama et al. 2006 for de- 
tails). The array is centered on the pixel containing the 
largest pulse height of an event. We combined the data 
from both modes as there did not appear to be any notice- 
able differences in the spectra obtained from the different 
modes. Data from version 0.7 of the pipeline reprocessing 
(Fujimoto et al. 2006; Mitsuda et al. 2006) were used for 
the analysis in the present paper. As is standard prac- 
tice, only ASCA grades 0, 2, 3, 4, and 6 were selected. 
Event files were cleaned to remove hot and flickering pix- 
els with the FTOOL cleansis as updated for the XIS. 
Data were excluded during satellite passages through the 
South Atlantic Anomaly (SAA), and for time intervals 
less than 256 s after passages through the SAA, using 
the T_SAA_HXD house-keeping parameter. Data were 
also rejected for Earth elevation angles (ELV) less than 
5°, Earth day-time elevation angles (DYE_ELV) less than 
20°, and values of the magnetic cut-off rigidity (COR) less 
than 6 GeV/c 2 . Version 0.7 of the reprocessing pipeline 
also includes CTI corrections based on available in-flight 
calibration and other refinements to the XIS energy scale. 
Residual uncertainties in the energy scale are on the or- 
der of 0.2% or less (or ~ 13 eV at 6.4 keV - see Koyama 
et al. 2006). We confirmed this from an analysis of the 
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Table 1. NGC 2992 Suzaku observation log. a Start and end times correspond to the time tags of the first and last photons 
respectively in the cleaned and filtered events files combined from XIS2 and XIS3. b Mean (0.5—10 keV) count rates and exposure 
times are for XIS2 and XIS3 per XIS. 



Obs 



Start 
(UT) 



End a 
(UT) 



Count Rate 6 
(ct/s/XIS) 



Exposure 
Time 6 (s) 



1 6/11/2005 14:16:51 7/11/2005 14:01:50 0.3888 ± 0.0027 34664.5 

2 19/11/2005 21:42:59 20/11/2005 23:24:45 0.5096 ± 0.0032 31641.5 

3 13/12/2005 10:15:37 14/12/2005 12:08:38 0.4408 ± 0.0026 41681.0 



onboard calibration line data taken during the NGC 2992 
observations (see section 2.3 for details). 

The cleaning and data selection resulted in net expo- 
sure times in the range 31.6 ks to 41.7 ks for the three ob- 
servations (see table 1). We examined XIS images made 
from the cleaned events files and found that the galaxy 
NGC 2993 (~ 3' from the core of NGC 2992) was clearly 
detected but its total count rate was less than 1% of the 
NGC 2992 count rate. Even though the half-power diam- 
eter of the XRT+XIS is - 2', NGC 2992 and NGC 2993 
are resolved. NGC 2993 was detected by Chandra with 
much better spatial resolution (see Colbert et al. 2005). 
We extracted XIS lightcurves and spectra using two sets 
of circular regions with radii of 2.6' and 4.35' centered on 
NGC 2992. The smaller region excludes NGC 2993, result- 
ing in negligible contamination of the NGC 2992 spectra, 
and the larger region corresponds to the size assumed by 
some of the standard XRT effective area files that have 
been released. The spectra from the smaller regions will 
be used for the principal analysis and those from the larger 
regions will be used to better estimate absolute fluxes (see 
below). We extracted background spectra using circular 
regions with radii of ~ 4' from source-free regions of the 
detectors. 

2.3. Instrument Response Files and the Mean XIS 
Spectrum 

For spectral analysis we used the response ma- 
trices (or "RMF" files) ae_xi0_20060213 .rmf , 
ae_xil_20060213c.rmf , ae_xi2_20060213 . rmf , and 
ae_xi.3_20060213.rmf for XIS0, XIS1, XIS2 and XIS3 
respectively (binned to the same number of channels 
as the spectra). The telescope effective area files (or 
"ARF" files) used were ae_xi0_hxdnom4_20060415 . arf , 
ae_xil_hxdnom4_20060415 . arf, 

ae_xi2_hxdnom4_20060415.arf and 
ae_xi3_hxdnom4_20060415.arf for XIS0, XIS1, XIS2, 
and XIS3 respectively. The ARF files are standard 
files released for the HXD-nominal pointing position 
and assume a circular extraction with a radius of 2.9'. 
Since our source extraction regions had a radius of 
2.6' we derived a mean correction factor by comparing 
broadband fluxes measured from the spectra extracted 
using 4.35' radii regions with those from our 2.6' radius 
region. The ARF files used for the spectra from the 
larger regions were ae_xiOJixdnom6_20060415 . arf , 
ae_xil_hxdnom6_20060415 . arf, 



ae_xi2_hxdnom6_20060415.arf and 
ae_xi3_b.xdnom6_20060415.arf for XIS0, XIS1, XIS2, 
and XIS3 respectively. Although the larger regions 
include NGC 2993, the contamination is small and the 
procedure is adequate for the present purpose. 

Some preliminary examination and spectral fitting re- 
vealed no significant evidence for variability in the total 
count rate or spectral shape within each of the three ob- 
servations. Although there was up to ~ 30% variability in 
net count rate between observations, we found no statisti- 
cally significant evidence for spectral variability between 
the observations (see table 1). We found that the four 
XIS sensors were consistent with each other to ~ 10% 
or better (within statistical errors) over the energy band 
0.5 — 10 keV. However, XIS0 showed a systematic dip in 
the ~ 5 — 7 keV band (i.e. in the region that Fe Ka line 
emission is expected) which was not apparent in the other 
three XIS sensors. Also, XIS1 showed a broad dip in the 
^2 — 5 keV band which was not apparent in the other 
three XIS sensors. Although both of these effects are at 
the level of ~ 10% or less, for the sake of maximum conser- 
vatism, in the remainder of this paper all of the spectral 
analysis will refer to a single XIS spectrum combined from 
observations 1, 2 and 3, using only XIS2 and XIS3 data. 
For analysis of the mean spectrum from XIS2+XIS3 (in- 
tegrated over the three observations), we combined the 
RMF and ARF files (four in total) into one response 
file with appropriate weightings, according to the expo- 
sure times of the data segments. The mean 0.5-10 keV 
background-subtracted count rate for this spectrum was 
0.4442 ± 0.0016 ct/s per XIS, and the net exposure time 
was 1.08 x 10 5 s per XIS. We restricted the analysis to the 
0.5-10 keV band in order to avoid background-subtraction 
systematics (in the 0.5-10 keV band the background con- 
stitutes ~ 3.4% of the total count rate). We obtained a 
flux correction factor for the size of the extraction region 
of 2.1% (as described above) for the mean spectrum. All 
fluxes and absolute line intensities quoted in this paper 
have been corrected (increased) by this factor. 

We measured the centroid energies of the Mn Ka 
lines from the calibration sources, using a background- 
subtracted spectrum combined from XIS2 and XIS3 data, 
summed over observations 1-3 (i.e. the same combination 
used for the analysis of NGC 2992), and summed over 
both calibration sources on each XIS. The expected ener- 
gies of the Mn Ka\ and Mn Ka-i lines are 5.89875 keV 
and 5.88765 keV respectively (Bcardcn 1967). Since the 
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Kol\ : Ka2 branching ratio is 2:1, the expected centroid 
energy is then 5.89505 keV. Using a single Gaussian we 
measured a centroid energy of 5.906 ± 0.003 keV (one- 
parameter, 90% confidence error). Therefore, the differ- 
ence between the measured and expected centroid energy 
is 11 ± 3 eV, consistent with the accuracy of the absolute 
XIS energy scale given by Koyama et al. (2006) as 0.2% 
(~ 13 eV) at the Fe Ka line energy. Both the residual 
CTI correction uncertainty and the spectral resolution are 
better at the centers of the detectors (i.e. in the region of 
the aimpoints for source observations) than at the detec- 
tor edges where the calibration sources are located so this 
is an over-estimate. We also measured the Mn Ka line 
width and obtained a FWHM of 161.5 eV (at an energy 
of ~ 5.9 keV). At 6.4 keV this corresponds to a FWHM 
of 168.4 eV (see Koyama et al. 2006), or - 7890 km s" 1 
FWHM. Since the FWHM resolution at 6.4 keV assumed 
in the response matrices is 133.8 eV (from direct inspec- 
tion), and measured line widths are a convolution of the 
assumed response with input lines, there is a residual 
FWHM width of y/(168A 2 - 133. 8 2 ) = 102.3 eV FWHM 
(or 4795 km s" 1 at 6.4 keV) that is not currently ac- 
counted for in the instrument response function. Again, 
the spectral resolution is better at the center of the detec- 
tors (by ~ 5 — 10 eV) than at the edges. 

2.4- XIS Contaminant Correction 

The XIS detectors have suffered a loss in low-energy effi- 
ciency due to absorption by a contaminant on the optical 
blocking filter. The degradation can be modeled using 
"effective" column densities of Carbon and Oxygen that 
can be calculated as a function of time since launch and 
position on each XIS, based on empirical fits to in-flight 
data obtained from non- varying astrophysical sources (see 
Koyama et al. 2006). We calculated the Carbon column 
density for each XIS and each observation of NGC 2992. 
For analysis of the mean spectrum described above we 
then calculated the average of the six column densities 
for XIS2 and XIS3, observations 1 to 3, and obtained 
3.4 x 10 18 cm -2 . Based on the current assumption that the 
contaminant is probably C24H38O4, the O column density 
was fixed at 1/6 of the value of the C column density. The 
C and O column densities were then combined with photo- 
electric absorption cross-sections from Balucinska-Church 
& McCammon (1992) since these were used to derive the 
empirical degradation formulae. The resulting contam- 
ination model, implemented using the varabs model in 
XSPEC, was included along with all astrophysical models 
of NGC 2992 (sec section 3). 

3. Spectral Fitting 

In the remainder of this paper we shall refer to anal- 
ysis of the combined spectrum from XIS2 and XIS3 av- 
eraged over all three observations (table 1). as described 
in section 2.3. We used XSPEC vll.3.2 for spectral fit- 
ting (Arnaud 1996), utilizing data in the 0.5-10 keV band. 
The instrument response files used were described in sec- 
tion 2.3. The spectrum was binned to ~ 29 eV per bin and 
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Fig. 1. (a) Ratio of the XIS2+XIS3 spectrum (integrated 
over observations 1 to 3) to a model consisting of Galactic 
absorption, an intrinsic power law continuum, and a uniform 
intrinsic absorber, (b) As (a), but now the model has addi- 
tional components corresponding to an unabsorbed power-law 
component (for example due to scattering), and narrow (unre- 
solved) Gaussian Fe Ka and Fe K/3 emission-line components. 
See section 3.1 for details. 

this resulted in > 25 ct /bin for the net spectrum over the 
0.5-9.7 keV range so we were able to use the x 2 statistic 
for finding the best-fitting model and for the statistical 
error analysis. Unless otherwise stated, all statistical er- 
rors given correspond to 90% confidence for one interest- 
ing parameter (A% 2 = 2.706). Galactic absorption with a 
column density of 5.26 x 10 20 cm -2 (Dickey & Lockman 
1990) was included in all of the models described here- 
after and its inclusion will be implicitly assumed. For the 
Galactic absorption, as well as any intrinsic absorption, 
we used photoelectric cross-sections given by Morrison & 
McCammon (1983). The correction for the XIS contami- 
nation degradation (as described in section 2.4) was also 
included in all of the spectral fitting described below. All 
astrophysical model parameter values will be given in the 
rest frame of NGC 2992. 

3.1. Preliminary Spectral Fitting and the Baseline Model 

Preliminary spectral fitting confirmed the complexity 
in the X-ray spectrum of NGC 2992 that has been estab- 
lished from observations with previous X-ray astronomy 
missions (see section 1). This complexity is illustrated in 
figure 1(a), which shows the ratio of the XIS data to a 
simple model consisting of a power-law intrinsic contin- 
uum and uniform intrinsic absorption fully covering the 
source. There are only three free parameters in this model, 
namely the normalization of the power-law continuum, its 
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photon index, T, and the column density, Nh- For this 
model x 2 = 1158.2 for 321 degrees of freedom. Relative 
to this simple model there is a large excess at soft X-ray 
energies below ~ 2 kcV and a very prominent Fe Ka emis- 
sion line. In addition, a weaker emission line is apparent 
at ~ 7 keV, which could be due to Fe xxvi Lya, Fe K/3, or 
a blend of both. From a detailed analysis described in sec- 
tion 5 we identify the weaker line with Fe K/3 and show 
quantitatively that contamination by any Fe xxvi Lya 
line emission is negligible. 

At least part of the soft excess in figure 1(a) is due to 
scattering of the intrinsic X-ray continuum into the linc-of- 
sight by an extended optically-thin zone in the active nu- 
cleus (e.g. Turner et al. 1998). This scattered component 
has been imaged in NGC 2992 using ROSAT and Chandra 
(see Colbert et al. 2005, and references therein) and in 
some other Seyfert 1.5-2 galaxies (e.g. Morse et al. 1995; 
Ogle et al. 2000; Young, Wilson, & Shopbell 2001; George 
et al. 2002; Iwasawa et al. 2003). The direct continuum 
along the linc-of-sight is obscured so even if only a small 
fraction of the intrinsic continuum is scattered, the soft 
X-ray spectrum can easily be dominated by the scattered 
component. Therefore, in the next step we added another 
power-law continuum component to the model, with the 
photon index tied to that of the direct continuum and a 
normalization that was allowed to float. We parameter- 
ized this normalization in terms of the fraction, f s , of the 
directly observed power-law continuum normalization. In 
the optically-thin limit the Thomson depth, r os , of the 
scattering region (assuming spherical symmetry) is sim- 
ply equal to f s . More realistically, f s = [ri/47r]r es , where 
[fi /Air] is the solid angle as a fraction of Air subtended by 
the scattering zone at the X-ray source. Only that part 
of the scattering zone visible to the observer should be 
included in this solid angle. In NGC 2992 an extended 
bi-conical outflow has been imaged (e.g. Allen et al.1999; 
Veilleux et al.2001), with opening angles of ~ 125° - 135°. 
Taking this angle to be 130°, we get [fl/Air] ~ 1 - cos 75° 
for the two cones together, or f s ~ 0.74r es . 

We also added two unresolved Gaussian components 
to model the emission lines in the Fe K band. The in- 
trinsic widths (<7jv) °f the two lines were tied together 
(in this case the two emission lines originate from the 
same atoms/ions). In the preliminary spectral fitting 
cn was fixed at a value (0.005 keV) much less than the 
XIS spectral resolution. The other parameters were the 
centroid energies of the two emission lines (E-p Ka and 
E^e K/3), and their intensities (I^e Ka and I-F e k/3 respec- 
tively). Thus, in this more complex model there were a 
total of eight free parameters (F, Njj, fs, Ep c k q , Ep c k/3, 
iFe Ka, Ipe k/3, and a continuum normalization). The 
best-fitting model gave % 2 = 400.6 for 316 degrees of free- 
dom. The ratio of the XIS data to this best-fitting model 
is shown in figure 1(b). The fit is still poor and statisti- 
cally significant residuals are still evident in the soft X-ray 
band below ~ 2 keV and in the Fe K band. Note that 
if wc included only one Gaussian component (to model 
the strongest emission line) we obtained x 2 = 418.8 for 
318 degrees of freedom, demonstrating that the weaker 



emission line is statistically significant (see section 5 for a 
more quantitative analysis). We then allowed the intrin- 
sic width (<tjv) of the Gaussian emission- line components 
to be free and re-fitted the model, obtaining another sig- 
nificant reduction in v 2 , which dropped to 375.4 for 315 
degrees of freedom. 

Motivated by the fact that an extended soft thermal 
emission component has been detected in NGC 2992 (see 
Colbert et al. 2005) and in other type 1.5-2 Seyfert 
galaxies (e.g. Ptak et al. 1997; Bianchi, Guainazzi, & 
Chiabcrge 2006, and references therein), we then added a 
component to model optically-thin thermal emission. We 
used the APEC model in XSPEC with the element abun- 
dances fixed at cosmic values (but later we allowed the Fe 
abundance to be free in order to derive constraints on it). 
The APEC model then involves two more free parameters, 
kT, and Lapec, the temperature and luminosity of the 
plasma respectively (although we used the normalization 
of the APEC model as the actual free parameter in the fits, 
not its luminosity). In all, there were now eleven free pa- 
rameters and we obtained \ 2 = 354.1 for 313 degrees of 
freedom. 

There still remained some residuals in the Fe K band in- 
dicative of additional, broader, line emission. This might 
be expected from an X-ray illuminated accretion disk 
whose inner radius extends to within ~ 10 or so grav- 
itational radii of the putative central black hole. Such 
relativistically broadened Fe K line emission has been 
observed in other AGN (e.g. Fabian & Miniutti 2005, 
and references therein). Indeed, it is possible that the 
Gaussian emission-line components in the above prelimi- 
nary model could be modeling part of the disk line emis- 
sion as well as the narrow core emission that might origi- 
nate from more distant matter. With CCD data the line 
core can be degenerate with part of the disk line compo- 
nent because both components can have significant emis- 
sion near the line rest-frame energy. Even though the 
Fe K line core may be the dominant component, it is in 
general necessary to simultaneously model both the broad 
and narrow components otherwise the inferred model pa- 
rameters for one component will be biased by the presence 
of the unmodeled component. For example, in the spec- 
tral fits described above with the Gaussian widths free, 
the best-fitting lines may be broader than the true width 
of the Fe K line core because they may be trying to model 
an underlying broad component as well. 

Therefore, we added Fe K line emission from a rela- 
tivistic disk around a (Schwarzschild) black hole using the 
diskline model in XSPEC (see Fabian et al. 1989). The 
diskline model has six parameters, namely the energy of 
the line in the disk frame, Eq; the inner and outer radii 
of the disk (i?; n and i? ut respectively, in units of gravita- 
tional radii, r g = GM/c 2 ); the power-law index, q, of the 
radial line emissivity law (r q ); the inclination angle of the 
disk normal relative to the observer's line-of-sight, 6> b s ; 
and the integrated intensity of the line, /disk- The overall 
Fe K line profile in the XIS data is already dominated by 
the line core and we found that the data for any resid- 
ual broad-line emission could not constrain all of the disk 



Precision Fe Ka and Fe K/3 Spectroscopy of NGC 2992 



7 



line model parameters. In particular, there is consider- 
able degeneracy between Eq and o bs- Therefore, we fixed 
Eq at 6.40 keV (corresponding to the rest-frame energy 
of Fe I Ka). In addition, we fixed R m and i? ut at 6r s 
(the radius of the last stable orbit around a Schwarzschild 
black hole) and 1000r g respectively. For self-consistency 
we also included a disk line component corresponding to 
Fe I K/3 with no additional free parameters: i?i n , i? ut, 
q, and 6* b s were forced to take the same values as the 
Fe Ka line, Eq was fixed at 7.056 keV, and the intensity 
of the Fe K/3 line was forced to be 13.5% of the Fe Ka 
line intensity (see section 5 for a detailed discussion on 
Fe K/3 line energies and the Fe K/3/Ka branching ratio). 
Although there may be no clear isolated detection of the 
broad Fe K/3 line, for ionization states of Fe I to Fe xvi 
it will be produced, and excluding it from the model could 
affect the values of the inferred parameters of the other 
Fe K line components. Fitting the XIS data with the Fe K 
disk line emission components included in the model re- 
sults in a drop in \ 2 of 19.7 (to 334.4 for 310 degrees of 
freedom) for an extra three free parameters (q, # bs, and 
^disk), indicating that a broad- line component is required 
by the data with a confidence level of more than 3cr. 

Finally, we added a component to model a Compton- 
reflection continuum using the model href 1 in XSPEC 
(see Dovciak, Karas, & Yaqoob 2004). Although only the 
^7 — 10 keV band XIS data arc able to constrain the re- 
flection continuum, if any of the Fe K emission line com- 
ponents originate in Compton-thick matter, a reflection 
continuum will be produced. Excluding a reflection contin- 
uum in the model could bias the best-fitting parameters 
of the Fe K emission lines. The reflection model involves 
only one extra free parameter, R, which is the normaliza- 
tion of the reflection continuum relative to the case of a 
steady-state X-ray continuum illuminating a slab subtend- 
ing a solid angle of 2ir at the X-ray source. The inclination 
angle between the disk normal and the observer's line-of- 
sight was tied to the disk inclination angle for the rela- 
tivistic Fe K emission-line components. Note that a disk 
geometry is likely to be inappropriate for any reflection 
continuum associated with Compton-thick distant-matter 
components of the Fe K line emission, but the data are 
not sensitive to differences in the reflection continuum due 
to geometry. When we fitted the above model to the data 
we found a best-fitting value of R = 0, and consequently 
no improvement in the fit. Nevertheless, we retained the 
reflection continuum in all of the spectral fitting analysis 
(with R free) , since it could potentially affect the statisti- 
cal error analysis for the other parameters. 

Thus, our baseline model consisted of fifteen free param- 
eters: two normalizations (one for the soft, optically-thin 
thermal continuum, another for the intrinsic power-law 
continuum), kT, V, N Hl f s , a Nl E Fc Ka , I Fc Ka , E Fc K( g, 
Ifq k/3, q, #obs, -fdiski and R. This model gave an excellent 
fit to the data (x 2 = 334.4 for 309 degrees of freedom), as 
can be seen from figure 2 (a) (which shows the best-fitting 
model overlaid on the XIS counts spectrum), figure 2 (b) 
(which shows the best-fitting model and some of the indi- 
vidual model components), and figure 2 (c) (which shows 
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Fig. 2. (a) The XIS2+XIS3 counts spectrum (integrated 
over observations 1 to 3) overlaid with the best-fitting baseline 
model (solid line). The inset shows a close-up of the region 
containing the Fe Ka and Fe K/3 emission lines. The C/O XIS 
contamination is included in the model, (b) The best-fitting 
baseline model (see section 3.1 and table 2). The soft X-ray 
optically-thin thermal continuum and the broad Fe Ka and 
Fe K/3 disk-line components are also shown separately. The 
C/O XIS contamination model has been removed in the plot 
but Galactic absorption is included, (c) The ratio of the data 
in (a) to the best-fitting baseline model. 



8 



T. Yaqoob et al. 



the ratio of the XIS data to the best-fitting model). The 
best- fitting parameters and 90%, one-parameter statisti- 
cal errors are shown in table 2. We discuss these results 
in detail in section 4. 



Table 2. XIS Spectral fitting results for NGC 2992. All val- 
ues and energy ranges in this table (except -Fo.5-2 kcV an d 
F 2 -w kcv) refer to the rest-frame of NGC 2992 (2 = 0.00771). 
Statistical errors are for 90%, one interesting parameter 
(Ax 2 = 2.706) and were derived with fifteen parameters free, 
with Tni Fo k« fixed at 0. Upper limits on the intensity and 
equivalent width (EW) of the Ni Kq line were derived with 
sixteen parameters free. 

a The upper limits on the widths of the line core of Fe Kq 
and Fe K/3 have not been corrected for the degradation of the 
XIS spectral resolution (see section 3.2). b Upper limits on 
R are for o b a fixed at the best-fitting value of 43° and for 
$obs = 60°, as indicated. c Observed-frame fluxes, corrected 
for the CCD contamination but not corrected for Galactic 
and intrinsic absorption. d Intrinsic, rest-frame luminosities, 
corrected for all absorption components. 



X 2 I degrees of freedom 
kT (keV) 

^APEC (10 40 erg s" 1 ) 

r 

N H (10 21 cm" 2 ) 

f s (scattered fraction) 

6 ohs (degrees) 

q 

-fdisk [Fe Ka] (10~ 5 photons cm -2 s -1 ) 
EW disk [Fe Ka] (eV) 

E N [Fe Ka] (keV) 

on (kcV) a 

FWHM (km s" 1 ) a 

7 N [Fe Ka] (1CT 5 photons cm~ 2 s" 1 ) 

EW N [Fc Ka] (cV) 

E N [Fe K/3] (keV) 

7 N [Fc K/3] (1CT 5 photons cm" 2 s" 1 ) 
EW N [Fe K/3] (eV) 

I N [Ni Ka] (1(T 5 photons cm" 2 s" 1 ) 
EW N [Ni Ka] (eV) 

Reflection fraction, b R(6 ohs = 43°) 
Reflection fraction, b R(6 ohs = 60°) 

Fo.5-2 kcv (10- 11 erg cm" 2 s" 1 ) c 
F 2 -io kcv (10- 11 erg cm" 2 s" 1 ) c 
£0.5-2 kcv (10 42 erg s- 1 )^ 

r nn42_-lVJ 



334.4/309 



h0.36 
-0.45 



1.18: 



1 ^gq+0-056 
1.0O9_q 027 

7 qq+0-56 

n7 o+0.021 
u - u ' °-0.021 



>31 

-1 5+ 5 ' 1 
-■-^-o.s 
1 q+o.5 
i - 9 -i.o 

h32 
-61 



118" 



6.407 



+0.007 
-0.007 

< 0.044 

<4850 



2.49 



+0.71 



163: 



-0.40 
,+47 
-26 

n 43+0.20 

u -^°-0.23 
27 +12 



^2-10 kcV 



(10 42 erg s- 1 ) c 



-14 

<0.29 

< 28 

<0.60 
<0.78 

0.16 
1.15 
0.70 
1.60 



3.2. Deconvolution of the Broad and Narrow Fe K 
Emission Lines 

Figure 3(a) shows a close-up of the Fe K region of the 



XIS spectrum of NGC 2992 in the form of the ratio of 
the data to a model consisting of a simple power-law con- 
tinuum with a fully-covering absorber (i.e. figure 3(a) is 
simply a "zoomed' region from figure 1(a)). Note that 
the Fe K/3 line is clearly seen here. A broad base to the 
prominent Fc Ka line core is also apparent. In section 3.1 
we found that when we modeled the underlying broad 
emission with a disk line, the Gaussian components that 
modeled the line core emission became narrower than for 
the case when the broad line emission was not modeled. 
In fact the best-fitting baseline model (table 2) indicates 
that the line cores of Fe Ka and Fe K/3 are unresolved by 
the XIS data. In order to better illustrate the Fe K line 
profile due to the putative rclativistic disk emission, we 
removed the disk line model components from the base- 
line model, froze the parameters of the thermal emission, 
reflection continuum, and Gaussian line core models and 
re-fitted the data. The resulting ratio of data to model 
is shown in figure 3(b) which illustrates the profile of the 
rclativistic disk component of the Fe K line. 

Next, starting from the baseline model again, with its 
fifteen free parameters (section 3.1) we constructed joint 
confidence contours of the equivalent width (EW) of the 
Fe Ka line core versus its intrinsic width (FWHM). The 
68%, 90%, 99% confidence contours are shown in fig- 
ure 3(c) which demonstrates that the 99%, two-parameter 
upper limit on the FWHM appears to be ~ 6300 km s _1 . 
However, as mentioned in section 2 the response matrices 
do not include the time-dependent degradation in the XIS 
energy resolution. In section 2.3 we showed that analy- 
sis of the Mn Ka calibration line data taken during the 
NGC 2992 observations gave a residual FWHM width of 
102.3 eV (or 4795 km s" 1 at 6.4 keV) that is not ac- 
counted for by the current response matrices. Since the 
observed line width from any observed source is a con- 
volution of the intrinsic line width and the instrumental 
broadening, the true 99% confidence, two-parameter up- 
per limit of the Fe Ka line core in NGC 2992 is then 
[6300 2 - 4795 2 ]^ - 4090 km s" 1 . Note that the narrow- 
line widths and their one-parameter statistical errors in 
table 2, the relativistic disk line parameters, and the con- 
tours in figure 3(c) have not been corrected for this effect. 
The one-parameter upper limit on the narrow-line width 
is similar to the residual error in the width of the response 
function (~ 102 eV FWHM or a ~ 43 eV). The disk line is 
measured to be a little broader than it really is so that our 
lower limit on the disk inclination angle is slightly over- 
estimated and the radial emissivity parameter, should 
be slightly smaller. 

We also constructed joint confidence contours of the in- 
tensity of the Fc Ka disk line versus the Fe Ka intensity 
of the narrow, distant-matter line. During this process, 
as for figure 3(c), all of the thirteen remaining parameters 
in the baseline model were free. The resulting contours 
are shown in figure 3(d), for 68%, 90%, 99%, and 3tr re- 
spectively. Remarkably, it can be seen that at a confidence 
level of greater than 3er, the disk and distant-matter Fe Ka 
lines are decoupled in the sense that the intensity of cither 
of the two line components is non-zero everywhere on the 
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Fig. 3. (a) A close-up of the XIS spectral region in figure 1(a) containing the Fe Ka and Fe K/3 lines. The vertical dotted lines 
correspond to the observed-frame energies of Fe I Ka and K/3 at 6.400 keV and 7.056 keV in the rest-frame respectively. All of 
the spectral analysis in the present paper refers to a mean spectrum constructed from XIS2+XIS3, integrated over observations 
1 to 3. (b) The ratio of the XIS spectral data to the baseline model re-fitted after removing the relativistic Fe Ka and Fe K/3 
disk-line components, whilst freezing the optically-thin thermal continuum, narrow Fe K line, and Compton-reflcction continuum 
parameters at their best-fitting values (see section 3.2 for details), (c) The 68%, 90%, 99% confidence contours of the equivalent 
width (EW) versus the intrinsic width (FWHM) of the narrow, distant-matter Fe Ka line based on the best-fitting baseline model 
(section 3.1). The contours have not been corrected for the degradation in XIS spectral resolution. After correction the 99% 
confidence, two-parameter upper limit is ~ 4090 km s _1 instead of ~ 6300 km s~ 1 (see section 3.2). (d) The 68%, 90%, 99%, and 
3a confidence contours of the Fe Ka relativistic disk-line intensity (/disk) versus the intensity of the narrow, distant-matter Fe Ka 
line (/narrow = Ijf) based on the best-fitting baseline model (section 3.1). It can be seen that even at a confidence level of 3cr, the 
disk and distant-matter Fe Ka lines are decoupled in the sense that the intensity of either of the two line components is non-zero 
everywhere on the 3<x contour. 
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3a contour. 

We investigated the possibility that part of the disk 
line may be an artifact due to the current XIS response 
matrices inadequately modeling the low-energy tail of the 
response function (see Koyama et al. 2006). Again we 
used the Mn Ka calibration line background-subtracted 
spectrum, integrated over observations 1-3, summed over 
XIS2 and XIS3, using both calibration sources on each 
XIS. We modeled the main peak with a disk line and 
Gaussian line, plus a simple power-law continuum. The 
disk- line model parameters (see section 3.1) i?i n , i? ou t 
were fixed at the values used for modeling NGC 2992 and 
the values of # bs, and q were fixed at the best-fitting val- 
ues shown in table 2. The disk line and Gaussian line 
intensities, as well as the peak energy and width of the 
Gaussian component were free parameters. By construct- 
ing confidence contours of the disk line intensity versus 
the Gaussian line intensity we established that the low- 
energy excess in the actual in-flight XIS response can only 
account for a disk line intensity that is 2% of the Gaussian 
line intensity, with a 99% confidence two-parameter up- 
per limit of ~ 6%. Yet for the actual NGC 2992 data we 
obtained a best-fitting Fe Ka disk line intensity that was 
76% of the line core (see table 2). Even the 3cr confidence 
contour of the disk line intensity versus the line core inten- 
sity (see figure 3(d)) only goes down to a value of ~ 8% for 
the ratio of the disk line to line core intensity. Therefore 
the low-energy excess in the tail of the XIS response does 
not have a significant impact on the measurements of the 
disk line in the NGC 2992 data or its decoupling from the 
Fe Ka line core. 

Decoupling of the Fe Ka components from the accre- 
tion disk and distant matter has not been possible with 
previous data for NGC 2992 (see section 4.4). It is even 
more remarkable that the decoupling has been possible 
with CCD resolution data and for a relativistic disk line 
that is not particularly impressive in terms of signal-to- 
noise ratio (figure 3(b)) or EW (see table 2). The reason 
can be traced to the fact that the narrow line has excel- 
lent signal-to-noise and a large EW. This in turn leaves 
less freedom for the disk line model. While it is possible 
that there may not be two distinct physical locations of 
Fe Ka line emission, a single line-emitting region would 
demand a peculiar line-emissivity law as a function of ra- 
dius. The strength of the line core compared to the broad 
base (which we model as a disk line) would demand that 
the cmissivity fall as a function of radius and then increase 
again. The width of the broad base of the Fe Ka line (see 
figure 3(b)) is - 33,000 km s" 1 FWHM, more than five 
times the upper limit on the width of the unresolved line 
core (figure 3(c)). Yet a uniform spherical distribution of 
matter gives line profiles which have much less contrast in 
the line width between the center and wings of the pro- 
file (e.g. see Yaqoob et al. 1993). The implications of 
the parameters derived from the dual Fe Ka line model- 
ing will be discussed further in section 4.4 and compared 
with published results on historical data. 



4. Discussion of Results 

In this section we discuss in detail the results that were 
obtained by fitting the baseline model (section 3.1) to the 
XIS data (see table 2). 

4-1- Soft, Optically-thin Thermal Emission 

We obtained a temperature and luminosity of the 
optically-thin thermal emission component of kT — 
0.656±g;glf kcV and - 1.2± 0.4x 10 40 erg s" 1 respectively. 
The temperature is typical of that found for a similar com- 
ponent in other type 1.5-2 AGN (e.g. see Ptak et al. 1997; 
Bianchi et al. 2006, and references therein). Note that 
the errors quoted here on the temperature are statistical 
only and do not include systematic errors, which may be 
larger. Chandra was able to image the soft X-ray emission 
better than Suzaku but was still not able to spatially sepa- 
rate it from other emission components due to contamina- 
tion from the PSF wings caused by the dominant nuclear 
emission (Colbert et al. 2005). Compounded with this 
were CCD pile-up issues, and lower signal-to-noise ratio 
of the data. Therefore Chandra obtained somewhat worse 
constraints on the temperature of the soft thermal emis- 
sion. Colbert et al. (2005) give kT = 0.51±°;?g keV and 
L = 1.4 x 10 40 erg s _1 and conclude that the bulk of the 
soft X-ray emission originates from within ~ 150 — 500 pc 
of the central engine. The temperature and luminos- 
ity of the extended soft emission component obtained by 
Chandra are both consistent with our respective values 
measured with Suzaku. 

Our Suzaku measurement of the temperature of the soft 
thermal emission in NGC 2992 is the most reliable to 
date. Elvis et al. (1990) reported detecting the soft X- 
ray emission component in NGC 2992 using the Einstein 
HRI, but it was only detected in one out of four quad- 
rants, the temperature was poorly constrained, and the 
luminosity was more than an order of magnitude greater 
than that obtained by Suzaku and Chandra. Although the 
soft X-ray emission imaged by the R OS A T HRI was spa- 
tially consistent with Chandra (see Colbert et al. 2005), 
the ROSAT HRI could not measure the temperature since 
it was not a spectroscopic instrument. The ROSATPSPG 
spectral resolution was too poor to measure the temper- 
ature. ASCA CCD spectra had better spectral resolu- 
tion but the signal-to-noise ratio of the NGC 2992 ASCA 
data is lower than that of the Suzaku data (see Weaver 
ct al. 1996). BeppoSAX could not constrain the soft ther- 
mal component because both the MECS and LECS in- 
struments had poorer spectral resolution than the Suzaku 
XIS and the useful MECS bandpass did not extend to 
below - 0.9 keV. 

4-2. Intrinsic Continuum and Absorption 

We obtained T = 1.5Tto.o3 f° r ^ ne power-law photon in- 
dex of the intrinsic continuum (see table 2). This is some- 
what less than that obtained from BeppoSAX (T ~ 1.7 
- Gilli et al. 2000) and Chandra (T ~ 1.8 - Colbert 
et al. 2005) observations. ASCA could not constrain T for 
the complex models required to explain the data (Weaver 
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ct al. 1996). Given the model-dependence of F and resid- 
ual instrument cross-calibration uncertainties of the order 
of Ar~ 0.1, we do not interpret the apparent differences in 
r amongst the different observations as necessarily having 
an astrophysical origin. 

We derived a column density of 8.0to'g x 10 21 cm -2 for 
the line-of-sight X-ray obscuration (see table 2) which cor- 
responds to a Thomson depth of ~ 0.006^q'qo4 (i.e. the ab- 
sorber is Compton-thin, at least in the line-of-sight). This 
column density is consistent with values obtained from 
previous observations with BeppoSAX (Gilli ct al. 2000) 
and ASCA (Weaver ct al. 1996). 

We tested the data to see if any additional Fe K ab- 
sorption was required by adding a simple absorption-edge 
model (zedge in XSPEC) with the threshold energy fixed 
at 7.11 keV (corresponding to neutral Fe) but allowing 
the optical depth at the threshold energy to be a free pa- 
rameter. All the other 15 free parameters in the model 
remained free. We obtained a best-fitting threshold op- 
tical depth of zero and a one-parameter, 90% confidence 
upper limit of 0.06. Since the optical depth at the Fe K 
edge threshold energy for the absorption already included 
in the model is ~ 0.008 (for N H = 8 x 10 21 cm" 2 and the 
Fe abundance used), we could only obtain a loose upper 
limit on the possible over-abundance of Fe of ~ 7.5. 

Table 2 gives the observed XIS fluxes (including the 
effect of Galactic and intrinsic absorption but not the 
CCD contaminant) in the 0.5-2 keV and 2-10 keV bands 
which were 0.16 and 1.15 x 10" 11 erg cm" 2 s _1 respec- 
tively. However, note that systematic uncertainties in 
the absolute flux calibration of X-ray detectors, as indi- 
cated by flux reference values quoted for the Crab Nebula 
observed by different instruments, are generally a few 
percent but can be as high as 10-20% in some cases 
(e.g. sec Jahoda et al. 2006, and references therein). 
Historically, the lowest and highest 2-10 keV observed 
fluxes for NGC 2992 have been measured by ASCA in 
1994 (~ 0.4 x lO" 11 erg cm" 2 s" 1 - Weaver et al. 1996), 
and by HEAO-1 in 1978 (~ 8.6 x 10~ n erg cm" 2 s" 1 - 
see Piccinotti ct al. 1982) respectively. This represents a 
factor ~ 21.5 in dynamic range. The Suzaku flux is at 
the lower end of the historically observed range, but it is 
nearly a factor ~ 3 higher than the historical minimum. 

Tabic 2 also gives the inferred XIS luminosities after 
correcting for absorption. These were calculated by set- 
ting the Galactic and intrinsic absorbing column den- 
sities to zero in the best-fitting baseline model (sec- 
tion 3.1). We obtained 0.7 and 1.6 x 10 42 erg s~ 4 in 
the 0.5-2 keV and 2-10 keV (rest-frame) bands respec- 
tively (assuming Hq = 70 km s~ 4 Mpc, A = 0.73, £1 = 1). 
The historical range in 2-10 keV luminosity is therefore 
~ 0.55- 11.8 x 10 42 erg s" 1 . 

4-3. Scattered Continuum 

Table 2 shows that we measured (7.3 ±2.1)% for the 
continuum scattered into the line-of-sight (as a percentage 
of the direct line-of-sight continuum before absorption). 
This translates directly into [n/4n]T cs = 0.073 ±0.021, 
where [Q/47r] is the fraction of the sky covered by the 



extended, warm, electron scattering zone, as seen from 
the X-ray source, including only that part of the scat- 
tering zone visible to the observer. Using [fi/47r] ~ 0.74 
for the bi-conical outflow in NGC 2992 (see §3.1), we get 
t os = 0.099 ± 0.028. The scattered fraction and optical 
depth cannot be easily compared with results from previ- 
ous observations. The ASCA data had very poor signal- 
to-noise and the model in Weaver et al. 1996 did not in- 
clude extended thermal emission. The BeppoSAX data 
did not have coverage of the soft X-ray band with suffi- 
cient resolution to measure f s (see Gilli et al. 2000). 

Colbert et al. (2005) determined from Chandra ACIS 
data that the scattering zone must be situated at more 
than ~ 150 pc from the nucleus since the spectrum from 
inside this region could be fitted with a power-law and 
absorption only. However, despite Chandra having much 
better spatial resolution than Suzaku, the Chandra data 
could not constrain the optical depth of the electron- 
scattering zone because the scattered continuum is spa- 
tially confused with the soft thermal emission and even 
with some of the direct central emission, due to the PSF 
wings. Spectroscopically, the Chandra ACIS data are in- 
ferior to the Suzaku data because they suffer from severe 
pile-up and have a much lower signal-to-noise ratio than 
the Suzaku data. For example, the 0.3-8 keV nuclear, 
and "near-nuclear" Chandra spectra analyzed by Colbert 
et al. (2005) had -6,744 and ~ 14,600 counts respectively, 
as opposed to a total of nearly — 10 5 counts (0.5-10 keV) 
in the XIS spectrum that we have analyzed here. Note 
also that pile-up in the Chandra data prohibited a reli- 
able measurement of the flux and luminosity of the nuclear 
source. Thus, even though the spatial resolution of Suzaku 
is not as good as Chandra, the high signal-to-noise ratio 
in the Suzaku spectrum and absence of pile-up problems 
allowed a better measurement of the scattered component 
and hence a better estimate of the Thomson depth of the 
scattering zone. 

4- 4- Fe K Emission Lines 

In section 3.2 and figure 3(d) we showed that the inten- 
sities of the narrow and broad Fe Ka emission-line com- 
ponents were decoupled. The best-fitting line parameters 
and statistical errors are shown in table 2. We obtained 
only a lower limit on the disk inclination angle, o b s > 31°. 
The actual value is somewhat dependent on the model 
used for the broad line (here it is the simple Schwarzschild 
black hole disk model which does not include light bend- 
ing - see Fabian et al. 1989). However, the overall roughly 
symmetric shape of the line (figure 3(b)) is a characteristic 
that excludes small inclination angles, independent of the 
detailed model, since smaller angles would give rise to an 
increasing asymmetry on the red side due to gravitational 
redshifting. We obtained an EW of H8TJ 2 eV for the 
broad line, which is consistent with that expected from a 
"cold" accretion disk illuminated by the observed direct 
X-ray continuum (e.g. see George & Fabian 1991). The 
index of the line emissivity radial power law, q = —1.5^08) 
is flatter than the average value for Seyfert 1 galaxies (e.g. 
Nandra et al. 1997a) but it is sensitive to the inner and 
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Fig. 4. (a) The 68%, 90%, and 99% confidence contours of 
^Fe K/3 versus Ep e K„(in the NGC 2992 rest- frame), mea- 
sured from the mean XIS spectrum (see section 4.4). The 
filled circles are theoretical rest-frame energies of the Fe K/3 
and Fe Ka emission lines from Palmeri ct al. (2003). (b) 
The 68%, 90%, and 99% confidence contours of Ip e k/3 ver- 
sus Tpe Ka(in the NGC 2992 rest-frame), measured from the 
mean XIS spectrum (see section 5). The lower and upper solid 
lines correspond to Fe K/3/Fc Ka branching ratios of 0.12 and 
0.135 respectively (see section 5). 

outer disk radius parameters, which we held fixed for the 
model-fitting (see section 3.1). We tried fixing q at the 
steeper value of —3 and obtained a lower limit on 8 Q ^ S 
and i?i n of 29° and 31r 9 respectively. 

For the centroid energy of the Fe Ka line core (which 
we modeled with a Gaussian component), we obtained 
6.407 ±0.007 keV (table 2). The statistical errors on this 
centroid energy are comparable to the residual systematic 
errors in the absolute energy scale of the XIS (on the order 
of 0.2% or ~ 13 eV at the Fe K line energy - see section 2.3, 
Koyama et al. 2006) . For the corresponding Fe K/3 line we 
obtained a centroid energy of 7.034^°;°^ keV. The 68%, 
90%, and 99% confidence contours of the narrow Fe K/3 
line energy versus the Fe Ka line energy are shown in fig- 
ure 4(a), and figure 4(b) shows corresponding contours for 
the line intensities of Fe Ka and Fe K/3. Note that both 
sets of contours were constructed from the baseline model 
(section 3.1) allowing, in each case, all thirteen of the re- 



maining parameters to be free. The measured energies and 
statistical errors on the Fe Ka and Fe K/3 line energies are 
a function of calibration uncertainties, Doppler and grav- 
itational shifts, and the ionization state of Fe. Further 
analysis utilizing information from the Fe K/3 line, includ- 
ing the Fe K/3/Ka ratio, is described in section 5. 

The measured intensity of the narrow (distant-matter) 
Fe Ka line, In = 2.5+° 'I x 10~ 5 photons cm~ 2 s _1 (see ta- 
ble 2), is a factor of ~ 20 times larger than the line inten- 
sity expected from a fully-covering Compton-thin spheri- 
cal shell with the cosmic abundance of Fe used in this pa- 
per, the measured line-of-sight column density (8.0+q 5 x 
10 21 cm -2 ), and the measured intrinsic 2-10 keV lumi- 
nosity (sec table 2). Using the expressions in Yaqoob 
et al. 2001 (see also Krolik & Kallman 1987), we obtained 
a predicted intensity of 1.1 x 10~ 6 photons cm~ 2 s _1 . 
Therefore, if the line emitter is Compton-thin, one or more 
of the following must be true: (i) the Fe abundance is 
higher than the assumed cosmic value, (ii) the column 
density out of the line-of-sight is greater than that mea- 
sured along the line-of-sight, and (iii) there are time delays 
between changes in the amplitude of the X-ray continuum 
luminosity and the response of the narrow Fe Ka line so 
that the observed line intensity may correspond to a time- 
averaged continuum illumination. On the other hand, the 
EW of 163^26 eV measured for the narrow Fe Ka line is 
also compatible with that expected from a line observed 
from the X-ray illuminated inner surface of a Compton- 
thick torus-like structure (although the historically aver- 
aged illuminating luminosity may need to be higher than 
the luminosity directly observed, e.g. see Makishima 1986; 
Ghisellini, Haardt, & Matt 1994). However, if the Fe Ka 
line core were produced in Compton-thick distant mat- 
ter we would expect to observe the so-called "Compton 
shoulder" on the red side of the line core. The Compton 
shoulder would have manifested itself as an asymmetry in 
the E-pe k/3 versus E-p e k q contours shown in figure 4(a), in 
the sense that they would be more elongated towards the 
lower energies (e.g. see Yaqoob et al. 2005 and references 
therein). Also, upper limits on the expected Compton- 
reflection continuum appear to be too low (section 4.5). 
Thus, the distant-matter Fe Ka line is likely to originate 
in Compton-thin matter. 

Table 2 and figure 3(c) show that the Fe Ka line core 
(narrow component) is unresolved by the XIS. We ob- 
tained a two-parameter, 99% confidence upper limit of 
4090 km s- 1 on the FWHM (or a N < 37 eV). These up- 
per limits on the line width have been corrected (as de- 
scribed in section 3.2) for the fact that the current instru- 
ment response matrices do not model the time-dependent 
degradation of the XIS spectral resolution (see Koyama 
et al. 2006). The measured upper limits on the line 
width can be compared with a theoretical estimate of 
FWHM ~ 760 yj (M$/r pc ) km s _1 (assuming a virial rela- 
tion and an r.m.s. velocity dispersion of \/3Vfwhm/2; e.g. 
see Netzer 1990). Here, Ms is the central black-hole mass 
in units of 10 8 M Q , and r pc is the distance (in parsecs), 
of the line-emitting structure from the central mass. The 
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mass of the central black hole in NGC 2992 has been esti- 
mated to be 5.2 x 10 7 M© from stellar velocity dispersions 
(Woo & Urry 2002, and references therein). Therefore, 
the measured upper limit on the FWHM corresponds to 
r pc > 0.018, or r > 21.3 light days. For comparison, broad 
Pa/3 line emission has been observed in NGC 2992 with 
a FWHM of - 1200 km s" 1 (Rix et al. 1990) indicating 
~ 250 light days for the location of the putative hidden 
BLR (see also Veilleux, Goodrich, & Hill 1997). 

It is pertinent to ask whether the narrow Fe Ka line 
core has varied in intensity compared to historical data 
(see Weaver et al. 1996; Gilli et al. 2000, and references 
therein). Given the large historical amplitude of variabil- 
ity of the continuum luminosity (greater than a factor of 
20), the response of the narrow Fe Ka line to the con- 
tinuum would give some clue about the size of the line- 
emitting structure. In order to answer the line variability 
question one needs to fit spectra with broad and narrow 
Fe Ka line models simultaneously, but this has not been 
done in the literature. We have fitted historical ASCA 
and BeppoSAX data with the dual Fe K line model used 
for the Suzaku data and have found that at the 99% confi- 
dence level, the ASCA (1994) and BeppoSAX (1997) data 
are consistent with no variability of the narrow Fe Ka 
line, and the BeppoSAX (1998) data cannot rule out vari- 
ability up to a factor of ~ 5. This is a factor of ~ 4 less 
than the continuum variability and supports an origin of 
the narrow Fe Ka line core in distant matter. At 99% 
confidence, the intensity of the broad (disk) Fe Ka line 
during the low-state ASCA and BeppoSAX (1997) obser- 
vations is consistent with no variability compared to the 
Suzaku measurement, but the high state BeppoSAX (1998) 
data cannot rule out a factor of up to ~ 15 variability. 
Therefore the relation between the broad line intensity 
and continuum luminosity remains open to future investi- 
gation. We note that in the 1997 BeppoSAX data (when 
NGC 2992 was in a near- minimum flux state), there were 
additional residuals at ~ 6.8 keV (rest-frame) relative to 
the dual-line model possibly due to emission from Fe XXV 
and/or Fe XXVI (see discussion in Gilli et al. 2000). 

4-. 5. Compton- Reflection Continuum 

If the broad Fe K line originates in an X-ray illumi- 
nated accretion disk then basic physics predicts a specific 
Compton-reflection continuum shape and amplitude com- 
mensurate with the Fe K line parameters (e.g. George & 
Fabian 1991; Magdziarz & Zdziarski 1995) and the disk 
ionization state. We used a simple disk reflection model 
(see section 3.1) to parameterize the total reflection con- 
tinuum from all possible locations (e.g. the accretion disk 
and any Compton-thick distant matter) because only a 
small portion of the XIS bandpass (^7—10 keV) can pos- 
sibly constrain it. The best-fitting value of the so-called 
reflection fraction, R, was 0. R is the normalization of 
the reflection continuum relative to that expected from a 
steady-state X-ray illuminated neutral disk subtending a 
solid angle of 2ir at the X-ray source. Although the XIS 
data cannot constrain the reflection continuum very well, 
we included it in the model because omitting it could po- 



13 



& 



O N - 

O 
ffl 

u 

c 
o 

o 

QJ 

K 




20 40 60 80 

e obs (degrees) 

Fig. 5. The 68%, 90%, 99% confidence contours of the rel- 
ative Compton-reflection fraction (R) versus the inclination 
angle (0 o b s ) of the disk normal to the observer's line-of-sight, 
based on the best-fitting baseline model (section 3.1). A value 
of R = 1 corresponds to a Compton-reflection continuum nor- 
malization expected from a disk subtending a solid angle of 2n 
at a central non-varying X-ray source. The inclination angle 
parameters of the disk for the Compton-reflection continuum 
and relativistic disk Fe Ka line model components were tied 
together. See section 4.5 for details. 

tentially affect the inferred values of some of the other 
model parameters, especially those of the broad disk line. 
We found that the upper limits on R were strongly depen- 
dent on the disk inclination angle, principally because of 
the angle-dependence of the relativistic line profile shape 
(the angular dependence of the reflection continuum shape 
in the restricted energy band is negligible). Therefore, we 
constructed joint 68%, 90%, and 99% confidence contours 
of R versus o bsj and these are shown in figure 5. It can be 
seen that the 90% confidence, two-parameter upper limit 
on R varies from ~ 0.6 to 0.9 as f? b s varies from its lower 
limit of 30° to ~ 80° respectively (the spectral fit is 
not well-behaved as we approach ~ 90°). For compari- 
son with values in the literature, the 90%, one-parameter 
upper limits on R are 0.60 and 0.78 for fixed inclination 
angles of o bs = 43° (the best-fitting value) and 6> b s = 60° 
respectively (see also table 2). 

We checked that the limits on the Compton-reflection 
continuum were not sensitive to background subtraction 
systcmatics by constructing a completely new background 
spectrum from rectangular regions from off-source loca- 
tions in XIS2 and XIS3 and summed over the three obser- 
vations. Note that the nominal background spectrum used 
thus far was constructed from circular regions centered at 
different positions and as far as possible relative to the 
new background regions (see section 2.2). With the new 
background spectrum in place, the baseline model with 
fifteen free parameters was re-fitted and there was no sta- 
tistically distinguishable change in any of the parameters. 
The one-parameter, 90% confidence upper limit on R for 
$obs = 43° remained at 0.60 and the corresponding upper 
limit for 9 ^ s = 60° decreased from 0.78 to 0.76. Note in 
figure 5 how the data strongly reject models with small 
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inclination angles (6* b s < 20° is ruled out at > 99% confi- 
dence). This is further confirmation that the narrow and 
broad Fc K lines arc decoupled because in poorer quality 
data a low-inclination angle disk line model can fit the 
narrow core of the Fe K line profile even if the line core 
has a non-disk origin. 

We also checked that the disk line intensity was not sen- 
sitive to the amplitude of the reflection continuum in the 
fit by constructing confidence contours of the disk line in- 
tensity, /disk, versus R (and reverting to the nominal back- 
ground spectrum for consistency with the earlier analysis). 
We found that that the 99% confidence lower limit on /disk 
was fairly flat and non-zero over the range < R < 3, and 
was never less than 0.4 x 10~ 5 photons enr 2 s _1 . 

Under the stated assumptions of our model we would 
expect R = 1. The one-parameter upper limits on R 
discussed above are not necessarily inconsistent with the 
measured EW (118±^ cV) of the broad Fc K line. The 
small bandpass leverage may bias R, and our model for the 
reflection continuum does not include relativistic smear- 
ing (which would tend to lower the effective value of R). 
The model also assumes a cold, neutral disk with a cos- 
mic Fe abundance. However, our upper limits on R do 
imply that there is unlikely to be a significant Compton- 
rcflcction continuum from the structure that produces the 
core of the Fe K fluorescent line emission (e.g. the putative 
obscuring torus) , which in turn implies that the structure 
is Compton-thin. 

Unfortunately the PDS data from both of the 
BeppoSAX observations, including the high-state obser- 
vation in 1998, were too noisy to obtain useful measure- 
ments or limits on a Compton-reflcction continuum (see 
Gilli et al. 2000). Although Colbert et al. (2005) claimed 
to measure a Compton-reflcction continuum in the off- 
nuclear region of NGC 2992 using Chandra, the data only 
extended up to 8 keV, the signal-to-noise was poor, and 
the data suffered from severe pile-up. The models fit- 
ted were complex in relation to the quality of the data 
and the uniqueness of the deconvolution of the different 
spectral components under such circumstances is an is- 
sue that was not addressed. In the future we expect the 
background modeling and subtraction systematics for the 
HXD onboard Suzaku to improve and we will revisit the 
challenge of measuring the Compton-reflection continuum 
in NGC 2992. Nevertheless, it is important to note that 
the detection of the relativistic disk line in the XIS data 
is insensitive to uncertainties in the reflection continuum: 
there was no preferred model solution in which a non-zero 
value of R forced the relativistic disk line intensity to zero. 

4-6. Upper Limits on the Ni/Fe Abundance 

We investigated Ni Ka line emission in the XIS data by 
adding another Gaussian line component to the baseline 
model (section 3.1). The intrinsic width of the additional 
Gaussian was tied to that of the Fe Ka and Fc K/3 lines. 
Since visual inspection revealed little evidence for Ni Ka 
line emission we fixed the center energy of the Gaussian 
at the expected (rest-frame) energy of 7.472 keV (Beardcn 
1967). Thus, only the intensity of the Ni Ka line was a 



free parameter, giving a total of sixteen free parameters 
in the model. The results are shown in table 2, and it 
can be seen that a statistically significant detection was 
not obtained. The decrease in the fit statistic was only 
A% 2 ~ 1.1. We obtained a 90% confidence, one-parameter 
upper limit on the EW of 28 eV. 

If the Fe Ka and Ni Ka lines are formed in Compton- 
thin matter, the ratio of their intensities depends princi- 
pally on the respective ratios of the fluorescence yields, 
the Kj3/Ka branching ratios, the photoelectric absorp- 
tion cross-sections, and the Fe and Ni abundances. Using 
Bambynek et al. (1972) for the fluorescence yields and 
K[3/Ka branching ratios, and Henke 1 for photoelec- 
tric absorption cross-sections, we evaluated equation (3) 
in Yaqoob et al. (2001) for Fe and Ni and obtained 
the theoretical ratio of line intensities: /Ni Kq//fc Ka = 
0.96(ANi/Ap e ). From the XIS data we constructed a joint 
90% confidence contour of /Ni kq versus /F e Ka and ob- 
tained an upper limit on /Ni Ka/lFo Ka of 0.176. Thus we 
obtained a 90% confidence upper limit on the Ni/Fe abun- 
dance ratio of 0.183. For comparison, the Ni/Fe abun- 
dance ratio from the Anders & Grevesse (1989) solar val- 
ues is 0.038 (i.e. a factor 4.8 smaller than our upper limit). 
If the Ka lines were observed in reflection from Compton- 
thick matter then the expected /Ni Kq//fc Kq ratio would 
be less than the Compton-thin case (due to absorption 
of the Ni Ka by Fe) and our upper limit on the Ni/Fe 
abundance would be larger. 

5. Precision X-ray Spectroscopy with the Fe Ka 
and Fe K/3 Lines 

In this section we will constrain the ionization state 
of Fc responsible for the narrow Fe K line emission in a 
robust manner by using as much information as possi- 
ble from the spectrum and by utilizing the latest atomic 
physics data and associated uncertainties. Firstly, we can 
make an immediate deduction from the very fact that the 
Fe K/3 fluorescence line is detected at all. That is, Fe 
must be less ionized than Fe XVII in the line-emitting mat- 
ter responsible for the narrow Fe K line emission because 
Fe XVII and ions with higher ionization states have no M 
shell electrons. We can constrain the ionization state fur- 
ther by taking advantage of the fact that we measured 
the Fe Ka and Fe K/3 line energies with sufficiently small 
statistical errors. We cannot do this directly from the 
measured line energies (e.g. using figure 4(a)) because 
they could be affected by residual XIS energy scale un- 
certainties, possible gravitational rcdshifts, and possible 
Dopplcr shifts. We will describe a method which enables 
us to place an upper limit on the dominant ionization 
state of Fe that is independent of the residual uncertainty 
in XIS gain, and does not require any knowledge of the 
possible gravitational and Doppler shifts. Our result will 
depend only on the uncertainty in the absolute offset in 
the XIS energy scale, on which wc can confidently place 
robust limits in general (Koyama et al. 2006), and for for 
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every observation using the onboard calibration sources. 
However, we first justify our confidence that the Fe K/3 
line suffers from negligible contamination from possible 
Fe XXVI Lya line emission. 

In figure 4(a) we showed the 68%, 90%, and 99% con- 
fidence contours of Ep e k/3 versus E-p e Ka measured from 
the XIS data. The fact that the contours are fairly sym- 
metric about the best-fit and are not elongated towards 
lower Fe K/3 energies implies that there is little or no 
blending from Fe xxvi Lya, whose rest-frame centroid 
energy would be at 6.966 kcV (e.g. Pike et al. 1996). We 
added another Gaussian line component to the baseline 
model (section 3.1) at a fixed energy of 6.966 keV, with 
its intrinsic width tied to that of the Fe Ka and Fe K/3 
lines and obtained a best-fit Fe xxvi Lya EW of and a 
90%, one-parameter upper limit of 18.5 eV. 

The 68%, 90%, and 99% confidence contours of Iy e k/3 
versus I-pe Ka measured from the XIS data are shown in 
figure 4(b). The lower and upper dashed lines correspond 
to ipc K,a/^Fo Ka ratios of 0.12 and 0.135 respectively. 
Palmeri et al. (2003) give a detailed discussion of the 
theoretical and experimental values of the Fe K/3/Fe Ka 
branching ratio. They compare theoretical values of the 
ratio in the literature for Fe I to Fe IX calculated using dif- 
ferent codes and assumptions and they compare these the- 
oretical values with several experimental measurements 
for Fe I. Palmeri et al. (2003) show that for Fe I, except 
for one experimental value, the remainder of the theoreti- 
cal and experimental values of the Fe K/3/Fe Ka branching 
ratio lie in the range ~ 0.12 — 0.135. Including the high- 
est experimental value, the range is ~ 0.12 — 0.145 corre- 
sponding to a dispersion of ~ 20%. For higher ionization 
states of Fe the branching ratio increases, and Palmeri 
et al. (2003) show that the theoretical value for Fe IX can 
be as high as ~ 0.17. The actual value of the observed 
ratio of the Fe K/3 to Fe Ka line depends a little on the 
optical depth of the medium in which the lines originate 
since the two lines have a slightly different albedo. 

New calculations for the energies of fluorescent Fe Ka 
and Fe K/3 transitions, fluorescent yields, and rates, for 
all charge states of Fe have been presented by Palmeri 
et al. (2003), Mendoza et al. (2004), Bautista et al. (2003), 
and Kallman et al. (2004). These works also include 
a comparison with previous calculations (and measure- 
ments, where appropriate) in the literature. In figure 4(a), 
for Fe I-Fe IX, we have plotted the expected Fe K/3 line 
energies against the expected Fe Ka line energies (solid, 
filled circles) , overlaid on the measured Fe K/3 line energy 
versus Fe Ka line energy confidence contours. The ener- 
gies were calculated from Palmeri et al. (2003), weighted 
according to the appropriate branching ratios for Ku\ and 
Ka>2 , and for Kj3\ and Kfo where necessary. As discussed 
in detail by Kallman et al. (2004) , and references therein, 
for ionization states higher than Fe ix the number of tran- 
sitions increases, and a single line centroid energy for com- 
parison with observations becomes less and less meaning- 
ful since the fluorescent line profile shape then depends on 
the details of the excitation mechanisms in play. One must 
then compute the detailed photoionization equilibrium in 



order to derive theoretical line profiles to compare with 
observations. From a comparison of results of theoretical 
calculations with laboratory measurements (solid state for 
Fe I, EBIT results for Fe ix), Kallman et al. (2004) assess 
the systematic error in the Fe K fluorescent line energies 
to be ~ 0.5 mA for Fe I and ~ 2 mA for Fe ix. We note that 
in eV, a systematic error of AA in mA, is AE(eV) ~ 3.3AA 
and A.E(eV) ~ 4.0AA at 6.4 keV and 7.0 keV respectively. 

It is important to note from the filled circles in fig- 
ure 4(a) that although the theoretical Fe K/3 line energy 
increases along the sequence Fe I to Fe ix, the Fe Ka line 
energy does not increase for these ionization states. This 
was pointed out by Palmeri et al. (2003) and Kallman 
et al. (2004) and the results are essentially a refinement 
of the results of Makishima (1986) which are universally 
quoted to describe a monotonic increase in the Fe Ka 
line energy with increasing charge state of Fe. However, 
Makishima (1986) never gave numerical values for the 
Fe Ka line energies, and for the ionization states Fe I 
to Fe ix it is not possible to extract the line energy val- 
ues from their plot with sufficient precision to compare 
with the values from Palmeri et al. (2003). Although the 
decrease in the Fe Ka line energy in going from Fe I to 
Fe ix is less than 10 eV, the XIS data for NGC 2992 are 
sensitive to these small energy differences. The statistical 
and systematic errors on the Fe Ka line energy measured 
from the NGC 2992 XIS data are comparable to the dif- 
ference in the theoretical values of the Fe I and Fe ix line 
energies. The values of the Fe Ka and Fe K/3 line energies, 
the associated statistical errors, and confidence contours 
measured for NGC 2992 are not only a function of the 
predominant charge state of Fe and residual calibration 
uncertainties, but are also affected by possible Dopplcr 
and gravitational shifts. In order to make inferences on 
the ionization state of Fe from the XIS data we must con- 
sider all of these factors, including uncertainties in the 
theoretical line energies, so we cannot simply deduce the 
ionization states allowed by the data from figure J^(a). 

In order to constrain the dominant ionization stage of Fe 
further, we can recast the Ep e k/3 versus E-p c Ka contours 
shown in figure 4(a). For a given ionization stage of Fe 
there are a pair of "true" (or theoretical) values of Ep c Ka 
and Ep c k/3- Suppose that the observed and true values 
of E-pc Ka and E-p c k/3 are linearly related by 

#i,Fe Ka (observed) = A + B x E itFc Ka (true) (1) 

and 

E iyFc K/3 (observed) = A + B x E uFc K /3 (true), (2) 

where i = to 25 represents the ionization state of Fe (i.e. 
the number of electrons missing in the ion relative to the 
neutral atom), A represents any residual calibration un- 
certainties in the offset of the absolute energy scale and B 
includes contributions from residual uncertainties in the 
XIS gain as well as possible gravitational and/or Dopplcr 
shifts. Certainly, over the restricted energy range being 
considered here, any higher order terms in the gain un- 
certainty are negligible (see Koyama et al. 2006), and the 
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Fig. 6. The 99% confidence contours of B versus A calcu- 
lated from the 99% contour of Ep e k/3 versus Ep e m fig- 
ure 4(a) (which was derived directly from the XIS spectrum of 
NGC 2992), for Fe I to Fe IX. The method of transformation of 
the Ep e k/3 versus Ep c k« contour is described in section 5. 
The Fe Kq and Fe K/3 fluorescent line energies were calcu- 
lated using atomic data in Palmeri et al. (2003) and are also 
shown as filled circles in figure 4(a). The variable A represents 
residual uncertainties in the offset of the absolute energy scale 
which is believed to be on the order of ~ 13 eV or less, but the 
dashed vertical lines are placed very conservatively at ±50 eV. 
B represents the product of any gravitational shift, Dopplcr 
shift, or residual uncorrected gain factors. Contours of B ver- 
sus A that do not have any portion lying inside A = ±50 eV 
are ruled out at 99% confidence. 

gravitational and Doppler terms only make a linear con- 
tribution. Therefore, the linear equations above are an 
excellent representation of the relations between the true 
rest-frame line energies and the final line energies mea- 
sured by the XIS. 

Since A and B are the same for the Fe Ka and Fe K/3 
lines, we can solve for A and B for a given ioniza- 
tion stage of Fe for any given pair of measurements of 
-E^Fe Ka (observed) and E^e k/3 (observed). Thus we 
can recast a contour of the measured energies Ev c k/3 and 
Eve Ka for a given confidence level into a contour of B ver- 
sus A for every i (i.e. we can generate up to 25 contours 
from figure 4(a)). Since we have independent information 
on the calibration of the absolute energy scale of the XIS 
data, we can then rule out any ionization states that yield 
B versus A contours that give solutions for A that are 
inconsistent with the known accuracy of the XIS energy 
scale. The absolute energy scale calibration is good to 
0.2% at the Fe Ka line energy and ~ ±5 eV below 1 keV 
(Koyama et al. 2006). 

The 99% confidence contours of B versus A calculated 
as described above are shown in figure 6 for Fe I to Fe IX, 
using the expected Fe Ka and Fe K/3 line energies calcu- 
lated from Palmeri et al. (2003). These energies were also 
shown as filled circles in figure 4(a). The vertical dashed 
lines in figure 6 represent the uncertainty in the absolute 
energy scale offset of the XIS (i.e. they enclose the possi- 
ble range of values of A). Although we have stated that 



this is nominally ~ ±13 cV or less, we will base our con- 
clusions on the very conservative assumption that \A\ is 
less than 50 eV. The vertical dashed lines are therefore at 
A = ±50 eV. It can be seen from figure 6 that the con- 
tours for Fe vii to Fe ix are ruled out regardless of the 
value of B. We can go further and take into account the 
uncertainties in the theoretical values of the Fe Ka and 
Fe K/3 line energies. The B versus A contours depend on 
Ep e K/3—E-pe Ka and we considered the two extreme cases, 
corresponding to the theoretical difference in these ener- 
gies hypothetically being in error by —2 mA and ±2 mA. 
We recalculated the B versus A confidence contours for 
each of these cases for Fe I to Fe ix. The cases for a the- 
oretical error of —2 mA and +2 mA gave contours which 
ruled out ionization states higher than Fe V and Fe vn re- 
spectively. Therefore, we can say that at 99% confidence 
the narrow core of the Fe Ka line in NGC 2992 originates 
in matter in which the predominant ionization state of 
Fe lies in the range Fe i-Fe VII. We emphasize that this 
conclusion requires no knowledge of the uncertainty in the 
XIS gain, or of any possible gravitational or Doppler shifts 
experienced by the photons leaving the source. 

6. Summary and Conclusions 

We have presented detailed X-ray spectroscopy in the 
0.5-10 keV band of the Seyfert 1.9 galaxy NGC 2992 with 
the Suzaku XIS. Here we summarize the results. 

1. The mean observed 2-10 keV flux of NGC 2992 dur- 
ing the three Suzaku observations in 2005 November 
and December was ~ 1.2 x 10~ n erg cm~ 2 s^ 1 , a 
factor of ~ 3 higher than the historic minimum flux, 
and a factor of ~ 6 lower than the historic maxi- 
mum. The absorption-corrected intrinsic 2-10 keV 
luminosity was 1.6 x 10 42 erg s _1 . 

2. We detected the extended, optically-thin soft X-ray 
thermal emission found with previous X-ray astron- 
omy missions, and obtained a temperature of kT = 
0.656^0 061 keV. This is the most precise measure- 
ment of kT for NGC 2992 to date although the errors 
are purely statistical and do not include any possible 
systematic effects. We obtained a luminosity for the 
soft emission component of~1.2± 0.4x 10 40 erg s _1 
which is consistent with that deduced from Chandra 
data (Colbert et al. 2005). 

3. We modeled the intrinsic X-ray continuum with a 
power law and obtained T = 1.57+o'o3 f° r ^ ne power- 
law photon index. 

4. We derived a column density of 8.0i"o'5 x 10 21 cm -2 
for the line-of-sight X-ray obscuration, which cor- 
responds to a Thomson depth of ~ 0.006^ ; 004 (i.e. 
the absorber is Compton-thin, at least in the line-of- 
sight). This column density is consistent with val- 
ues obtained from previous BeppoSAX and ASCA 
observations. 

5. We estimated r cs = 0.10 ±0.03 for the optical depth 
of the extended, warm, electron scattering zone 
which is responsible for scattering (7.3 ±2.1)% of 
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the intrinsic hard X-ray continuum into the line-of- 
sight. 

6. We detected broad and narrow Fc K line emis- 
sion. When we modeled the broad and narrow lines 
with line emission from a relativistic disk around 
a Schwarzschild black hole and a Gaussian profile 
respectively, we found that the intensities of the 
broad and narrow lines were decoupled at a con- 
fidence level of greater than 3a. We obtained an 
EW of 118^61 cV, a radial line-emissivity law of 
r -i.s(-o.8,+s.i) ! and a lower limit of 31 o Qn thc disk 

inclination angle from the disk-line model. The ac- 
tual values are model-dependent, however. 

7. Thc narrow Fe K line core EW was 163+26 eV - The 
line core was unresolved by the XIS and wc obtained 
a 99% confidence, two-parameter upper limit on the 
FWHM of 4090 km s^ 1 , after taking into account 
the time-dependent degradation in the XIS spec- 
tral resolution (as discussed in section 2.3 and sec- 
tion 3.2). 

8. Thc absolute flux in the narrow line is a factor of 
more than 20 larger than that expected from a fully- 
covering spherical shell illuminated by the measured 
intrinsic 2-10 keV luminosity, having a cosmic abun- 
dance of Fe and a uniform column density equal 
to that measured in the line-of-sight. Therefore, 
roughly speaking, the product of the solid angle sub- 
tended by the line-emitting structure at the source, 
its angle- aver aged column density, its Fe abundance, 
and the mean historical 2-10 keV luminosity aver- 
aged over timescales of the order of the structure 
light-crossing time, must be > 20. 

9. We detected the Fe K/3 fluorescent line correspond- 
ing to the Fe Ka line core. Taking into account 
residual uncertainties in the XIS energy scale cali- 
bration and conservative uncertainties on the theo- 
retical energies of the Fe Ka and Fe K/3 transitions 
for different ionic species of Fe, we were able to con- 
strain the predominant ionization state of Fe respon- 
sible for producing the narrow Fe K line core to be 
Fe i-Fe vn. Ionization states of Fe vin or higher 
are ruled out at 99% confidence. This is the most 
precise determination to date of the ionization state 
of the matter emitting the Fe K line core and was 
achieved using a method that requires no knowledge 
of any uncertainties in the XIS gain and any possible 
gravitational and Doppler shifts. 

10. We were able to decouple the disk and distant- 
matter Fe Ka emission lines because the continuum 
flux was low, giving a large EW for the line core, 
with good signal-to-noise ratio. The large core EW 
also enabled the Fe K/3 line to be measured with 
a good signal-to-noise ratio, giving redundant in- 
formation on the width of the distant-matter lines. 
This demonstrates the importance of measuring the 
distant-matter Fe K line emission parameters in or- 
der to measure those of the relativistic Fe K line. 

11. Although the XIS data do not have the bandpass 
to measure the Compton-rcflcction continuum as- 



sociated with the broad Fe K line, we included it 
in all of the spectral analysis with a parameterized 
normalization, R, relative to that expected from a 
steady-state illumination of a disk subtending a solid 
angle of 2ir at the X-ray source. We obtained upper 
limits on R as a function of thc disk inclination angle 
(constrained by the relativistic Fe K line). We found 
(at 90% confidence for two parameters), R < 0.9 for 
an inclination angle of 80°, going down to R < 0.6 
when the inclination angle is at its lower limit of 
31°. If all of this reflection is associated with the 
broad Fe K line, which presumably originates in a 
Compton-thick disk, there is no room left for reflec- 
tion from distant matter. This implies that the nar- 
row Fe K line core is likely to originate in Compton- 
thin matter, supported by the lack of a Compton 
shoulder on thc Fc K line core. Although thc EW of 
the Fe K line core is too large for steady-state illu- 
mination of Compton-thin matter, several possible 
factors can account for it (see [8] above). 

In the future we will refine our modeling by utilizing 
the Suzaku HXD data when the background model has 
achieved the goal of 1% systematic uncertainties. The 
HXD data will better constrain the Compton-reflection 
continuum components. We will also improve model con- 
straints by including data from XIS0 and XIS1 when 
the calibration has advanced further. Future observa- 
tions with higher spectral resolution, such as gratings or 
calorimeters, will be able to resolve the narrow core of 
the Fe Ka line and thereby better constrain the location 
of the line-emitter. In addition, improved throughput will 
enable better measurement of the profile of the broad Fe K 
line. Further monitoring of NGC 2992 to measure vari- 
ability in thc Fc K line components and in the reflection 
continuum as the continuum undergoes large excursions 
in amplitude promise to reveal important clues on thc ac- 
cretion process and on the structure of the central engine 
beyond the accretion disk. 
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